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Abstract 

We study in detail the contribution of heavy composite Majorana neutrinos 
to neutrino-less double beta decay (On/3/3). Our analysis confirms the result 
of a previous estimate by two of the authors. Excited neutrinos couple to 
the electroweak gauge bosons through a magnetic type effective Lagrangian. 
The relevant nuclear matrix element is related to matrix elements available 
in the literature and current bounds on the half-life of On/3/3 are converted 
into bounds on the compositeness scale and/or the heavy neutrino mass. Our 
bounds are of the same order of magnitude as those available from accelerator 
experiments. 
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I. INTRODUCTION 


Neutrinoless double beta decay (0z//3/3), see Fig. 1, is certainly one of the more interesting 
non-accelerator processes that are presently being searched for. The interest in this process 
stems from the fact that its observation would undoubtedly signal lepton number violation, 
and at the same time would shed light onto the nature of the neutrino, one of the most 
elusive elementary particles. For these reasons it has received considerable attention both 
from the nuclear and the particle physics community |^. 

In the standard model, QuPP can only proceed if the neutrino is of Majorana type and 
has a non zero mass. A number of mechanisms studied in models beyond the standard 
electroweak theory |^, have verihed that neutrino-less double beta decay is a very good 
probe of physics beyond the standard model. The experimental lower bound on the lifetime 
of the decay has been used to obtain constraints on the scale of new physics. 

Recent work along these lines include: (i) an investigation of new super-symmetric con¬ 
tributions from R-parity violating MSSM [Q shows that constraints on parameters of the 
model from non-observation of 0z//3/3 are stronger than those available from accelerator ex¬ 
periments; {ii) a detailed analysis of the contribution from left right symmetric models [Q; 
{ill) a study of the effective low energy charged current lepton quark interactions due to the 
exchange of heavy leptoquarks |^. The phenomenology of Majoron models has also been 
studied in detail @]. 

Panella and Srivastava 0] were the hrst to show that the compositeness scenario can give 
an additional contribution to 0z//3/3 and derived bounds on the compositeness parameters 
from the non observation of Ouf]/]. They explored phenomenologically the idea that the 
excited state of an ordinary neutrino might be a heavy Majorana neutral with a mass 
ranging from a few hundred GeV up to a TeV. However, the nuclear aspect of the calculation 
was treated only approximatively: time ordering of the hadronic charged current was ne¬ 
glected and an upper bound for the nuclear matrix element was used in deriving constraints 
on the compositeness parameters. Use of the HEIDELBERG-MOSCOW (3(5 experiment 
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lower bound |§] on the half-life of the decay ™Ge —Se -1- 2e“, yielded the following con¬ 
straint [| on the scale (Ac), the heavy neutrino mass (M^v), and the dimensionless coupling 
constant / |^: 


I/I <3.9 


Ac 

ITeV 


ITeV; 


( 1 ) 


Apart from the obvious desire to improve on the above mentioned approximations the 
main motivation for the present work is twofold. Firstly, after the completion of the work 
of ref. 0], appeared a related work by Takasugi [§] who considered the same problem 
via the exchange of an heavy composite Majorana neutral) but arrived at quite different 
conclusions P [|: 


o Ac /lTeV\V2 

< 5.8 X —— . 

1 TeV V Mjsf J 


( 2 ) 


In view of this discrepancy it is, of course, mandatory to investigate further the problem in 
order to understand if reliable constraints on the compositeness scale from are given 

by Eq. (|I|) or by Eq. (1). 

Secondly, there is a recent claim by the CDF collaboration of a possible signal of com¬ 
positeness in high energy proton-proton collisions |^: the measured cross section for large 
transverse energy jets is signihcantly higher than predictions based on perturbative QCD 
calculations to order O(a^); a compositeness scale of Ac = 1.6 TeV is suggested by the CDF 
study. Were this claim to withstand further data and analysis (such as angular distribution 
of dijets presentely underway), the interest in new physics effects arising from a compos¬ 
ite scenario will undoubtly increase enormously. If so, low-energy processes, such as Oz//?/?, 
could play a complementary role and hence are worth further investigation. 

In this paper we present a detailed analysis of the composite Majorana neutrino contri¬ 
bution to 0z//d/9; we show that: {i) the peculiarity of the dimension Eve effective coupling 


^the numerical value used in ref. Q was: > 5.1 x lO^^yr. 

^ the numerical value used in ref. Q was: > 5.6 x lO^^yr. 
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{o'liu) shows up in giving a larger than usual importance to the high-energy behaviour of the 
hadronic current correlation function; [ii) the nuclear matrix element is calculated exactly 
since it can be related to matrix elements already known; [Hi) the results of Panella and 
Srivastava [0 remain essentially unchanged; (iv) the calculation by Takasugi presented in 
ref. is not consistent, and its conclusions are not correct. 


II. EFFECTIVE LAGRANGIANS FOR COMPOSITENESS 


The idea that at an energy scale Ac, quarks and leptons might show an internal struc¬ 
ture has been around for quite some time 0. Although many models describing quarks 
and leptons in terms of bound states of yet more fundamental entities (preons) have been 
proposed, so far, no consistent dynamical composite theory has been found. Phenomenolog¬ 
ically however this idea can be probed by observing that one natural, model independent, 
consequence of compositeness is the existence of excited states of the ordinary fermions with 
masses at least of the order of the compositeness scale. 

We review here, for the reader’s convenience, the effective interactions used in the litera¬ 
ture to describe possible manifestations of lepton and quark substructure. A more detailed 
discussion of compositeness phenomenology can be found in refs. |0,0. 

Effective couplings between the heavy and light leptons and quarks have been proposed, 
using weak isospin {Iw) and hyper-charge (V) conservation. Within this model, it is as¬ 
sumed that the lightness of the ordinary leptons could be related to some global unbroken 
chiral symmetry which would produce massless bound states of preons in the absence of 
weak perturbations due to SU{2) x U{1) gauge and Higgs interactions. The large mass of 
the excited leptons arises from the unknown underlying dynamics and not from the Higgs 
mechanism. 

Assuming that such states are grouped in SU{2) x U{1) multiplets, since light fermions 
have Iw = 0,1/2 and electroweak gauge bosons have Iw = 0,1, only multiplets with Iw < 
3/2 can be excited in the lowest order perturbation theory. Also, since none of the gauge 
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fields carry hyper-charge, a given excited multiplet can couple only to a light multiplet with 
the same Y. 

In addition, conservation of the electro-magnetic current forces the transition coupling 
of heavy-to-light fermions to be of the magnetic moment type respect to any electroweak 
gauge bosons |]^. In fact, a 7 ^ transition coupling between e and e* mediated by the 
and gauge fields, would result in an electro-magnetic current of the type ~ ' 0 e* 7 ^t/’e 
which would not be conserved due to the different masses of excited and ordinary fermions, 
(actually it is expected that rrie* 3> me). 

Let us here restrict to the hrst family and consider spin-1/2 excited states grouped in 
multiplets with Iw = 1/2 and Y = —1 (the so called homodoublet model |T^), 


L = 




which can couple to the light left-handed multiplet 


= 


.ei, 


1 - 75 / 

2 U, 


( 3 ) 


( 4 ) 


through the gauge helds andB^. The relevant interaction can be written |T^ in terms 
of two new independent coupling constants / and /': 


iVp Z 

, g'f 


Ar 


-La^jL]-d'^B^ + h.c. 


( 5 ) 


where r are the Pauli SU (2) matrices, g and g' are the usual SU{2) and U{1) gauge coupling 
constants, and the factor of — 1/2 in the second term is the hyper-charge of the U{1) current. 
This effective Lagrangian has been widely used in the literature to predict production cross 
sections and decay rates of the excited particles ||T3|JT5| , pT . 


The extension to quarks and strong interactions as well as to other multiplets and a 
detailed discussion of the spectroscopy of the excited particles can be found in the litera¬ 
ture HI? . 


Here, let us write down explicitly the interaction Lagrangian describing the coupling of 
the heavy excited neutrino with the light electron, which will be slightly generalized in the 
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following section in order to discuss bounds on the compositeness effective couplings from 
low-energy, nuclear, double-beta decay: 


^eff - 




. 1-75 


V 2 A, 


( 6 ) 


III. NEUTRINO-LESS DOUBLE BETA DECAY (Ou00). 


The transition amplitude of Oiy00 decay is calculated according to the interaction La- 
grangian: 


r — 




V2 


+ cos9c\ x) + h.c. 


( 7 ) 


where we have generalized the interaction in Eq. (P) to include right-handed couplings (in 
order to allow comparison with other models than the homo-doublet one), although we will 
assume chirality conservation i.e. {rjL^rjR) = (1, 0) or (0,1); R= (1-|- 75 ) 72 , L = (1-|- 75 ) 72 , 
9c is the Cabibbo angle (cos^c = 0.974 ) and is the hadronic weak charged current. 

We have: 


2 / 9 \4 f f 




‘^‘1 J 4 . J 4 


[27ry 


crxcTy exp[—iq ■ {x — y)] x 


4 X ^(1 - Pi2)i’{p2)(Tf,x{VLL + prR) {plL + yRR)aup'ip (pi) x 

exp[i(pi ■x + p2-y)]{q- Pi)^(g + ^2)^1 


+ M]y iC( 

- - ^{plL + PRR)a,4' 

q 

{F\T[jj:{x) javw) 


( 8 ) 


[(g - Piy - M^] [{q + P 2 y - M^] 
where (1 — Pi 2 )/V^ is the antisymmetric operator due to the production of two identical 
fermions. Now let us change variables of integration: 


X = z H- 2 


with d^x d'^y = d^u. 


( 9 ) 


2 / = ^-f 

In addition we make the ansatz that the hadronic current be given by the corresponding 
sum of the nucleonic charged current: 
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i 

This implies: 


{F\TlJj:(x) rM]\I) = exv\i(PF - Pi) ■ »)] (F|rK(i - y) J^(0)]|/) 


We have: 


• pi-x + p2-y = {P1+P2) ■z + {pi- P2) ■ u /2 

• {plL + PrR){4 + Mn){plL + PrR) = (pIL + pI^R)Mn if PrPl = 0. 

. {g/{2V2)r = M^GI/2 


Thus we arrive at: 


S„ = {cosec)^Gl U-^ I 


dS 


d'^zd'^u exp{—iq ■ u) x 


Mn ^(1 - Pi 2 )i}{p 2 )a^xCTup{vlL + plR)ij^{pi) x 
exp[i z- {P1+P2 + PF- Pi)] exp [i (m/ 2 ) ■ (pi - P2 - Pf + Pi)] x 
W Pi) W + PP 1,^ _ ) 

The integration over d'^z gives the energy-momentum conservation and if we define: 


Sfi = ( 27 r)^( 5 ^(p 7 -PF-Pi- P2) r 


fi 


^eff — COS 9 c Gf (//Ac) 


we obtain: 


Tfi = (1 - P12) J d\ exp[-i{q - pi) ■ u] x 

Mjv %i}{p2)a^x(yvp{rilL + plR)ij^{pi) {q - pi)^{q + ps)^ x 

_ (f|t[j/‘(m) jno)]|/) _ 

[(g _ p^)2 _ M^] [(g + P2)2 - M^] (g2 _ M^) 


Next, we neglect pi and p 2 everywhere with respect to q except in the electronic 
functions. Using the identity: 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 

wave 
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(1 - Pia) ij{p2)af,xaup{vlL + r]\R)i)^{pi) = '4){p2){(J^,\, (^upjivlL + p\R)ip^{pi) 

= 2^{p2){ppuVxp - VwVxu + i-/5e^xvp){plL + p\R)i)^[pi] (14) 


we find: 


2G^ 

Tf. = -^iip^)(r,lL + r,lR)ifi(t,i)M^M^ 


d*q (j'rr - qV) Wv(«) 


(q^ - + k)‘^{q^ - M% + ie) 

( 16 ) 


where we have dehned: 


= / d'x exp(-tq ■ x) {F\T[r^(x) 4"(0)]|/). 


( 16 ) 


Eq. 


gives: 


Tf^ = 


2Gl 


eff 


i^{P2){vlL + PRR)^lJ^{pl) MnM^ 


(Rq 


X 


f dqo -q^q^[Wqj - p^.iWoo + W^)] + ql 


(17) 


. 2tt {ql -uj% + ie) {ql - o;^ + ie)^ 

where we have dehned: ujw,n = and terms in go^i have been dropped because 

they do not contribute to O’*" —O'*' transitions. 

Inserting a complete set of intermediate states one can cast Wp,qy in the form: 


Wpu{q) = {-i) / d^x exp{iq ■ x) ^ 


(F|J/((x)|X)(X|J,^(0)|/) , {F\J!:{0)\X){X\j;i{x)\I)^ 


qo - Ep + Ex - ie 


+ 


-go - Ej + Ex - ie 


i 

( 18 ) 


Using some known results on matrix elements of one particle operators, the quantity 
IU^i,(g) can be readily evaluated. Note that the sum over intermediate states includes an 
integration over the continuous part of the spectrum, namely the center of mass momentum 
P: J2x {27 i)~^ f d?P Yxn- This integration can be carried out analytically. 

The energy of a state |X) is Ex = Ecm{P) + Gi where Ecm{P) is the energy of 
the center of mass translational motion and e^, is the excitation energy [jT^. {Ecm{P) = 
\JP"^ + (rUpA)^, with rup the proton mass and A the mass number of the nucleus). 













In the closure approximation (routinely applied in double beta calculations), the excita¬ 
tion energy e^, of the intermediate state is replaced by an average value and the sum over 
the discrete part of the intermediate states is performed. 


The center of mass motion (R = (1/A) ) can be separated out so that we have [IK 


i, = ri-R = 0 

i 

In terms of the new coordinates (^J we have: 

i’iXi, .,r^) =exp(sP-fl) <*Ki,. 


(19) 


( 20 ) 


Using this for the calculation of one body operator matrix elements [|^ one obtains (for 
states of dehnite center of mass momentum): 


{F|0(r)|.Y) = y {{FI exp|i (Py - P,) ■ (r - ^,)l 0®(Px - Pf) \X}} 


( 21 ) 


where we denote with ((...)) the matrix element in the space of the A —1 relative coordinates. 
After integrating over the center of mass motion, we obtain: 

2iA 


= 2 ■ ^ki)jl!"\-q)jj!;\q)\i)) 

qo - n, -\-ze 

(in the limit pp ^ pj ^ 0). We have also used the approximation ep k, ej : 


( 22 ) 


^ — Cn — ^7 + EcM{q) — Ecm{^) 


(23) 


Since on the average |g| ~ 40 MeV, we may conclude that the center of mass motion does 
not give any appreciable contribution to A: 


A Ri Cn — 6/ 10 MeV 


(24) 


Next, we use the non-relativistic limit for the nuclear current [19 

= /a(9t X 


(fc) 

gvrX’ 


if p = 0 
-gAT+^^\cFk)i if p = 7 
where a*, is the spin matrix of the fc-th nucleon, and a nucleon form factor. 


( 26 ) 
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(26) 



with niA = 0.85 GeV, has been introduced to account for the hnite size of the nucleon. The 
latter is known to give a sizable contribution for the heavy neutrino case. 

Within the above approximations, we are thus led to: 



9 \ {{F\ Y, exp(^g ■ ru) rf Vf (afc)i(cTi)i|/)) (27) 


These expressions for W^i, are standard and the closure approximation has been routinely 
applied in double beta calculations p. 

IV. DECAY FORMULAE AND NUCLEAR MATRIX ELEMENT 

We can now proceed to calculate the Oz//?/? decay amplitude using the results of the 
previous section in Eq. (pTf) . 

Dehning: 





1 


+ F) {qi -u;w + F) 



(28) 


we can write: 



kl 


-sW-) ■ £'(0 r(92)||/» X fl(q^) 


(29) 


The integrals J and J' are: 
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( 30 ) 


1 A + UJy\r + Ct^TV 

2 U}\yU}n {m^W + ^n') (A + U}-\y') (A + l^Tv) 


2 (ujw + ^n) (A + ujw) (A + ujn) 


We note that between the three scales A, ujy/ and a ;at involved in the problem we have 
the following ordering: 


A U)wr Ci^jv 


(31) 


Using Eq. 


one Ends: 


I 

I' 


Atuvi/ 

, ,4 , ,2 


W 


l + |A + ofA>= 

2 Ctivv V 


(32) 


Using these results into Eq. (|^) we hnd: 


- 

N 


HP2){vIL + ri\R)i)^{pi) g\ X 


^/2M, 


{{F\Y. '^+^^+{Dij{rki) 


kl 


(fc) (0 (9v 

^9 A 


-B(r*,)<7<‘l.cr<"||/)) 


(33) 


where Dij{rki) and D{rki) are given as: 


Dijir-ki) = [ 7 ^-^ exp {iq ■ r^/) 


(27r) 


(1 + qym\y 


(Pq .. , AiU^y 


/ d q 

exp{iq-rki) 


(1 + qym\Y • 


(34) 


Explicitly : 




m\ Rq 


A-kRq Vki 


SijFAixA) - 


{rki)i {rki)j 


Fb{xa) 


kl 


AMwm\ Ro ^ , 

D{rki) = —r-B - Fn{xa) 

AtvRo Tkl 

where Rq = tqA^R (ro = 1.1/m) , xa = rnATki and : 


( 35 ) 
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Fa{x) = — exp(-x) + x) 


Fb{x) = — exp(-x) 


Fn{x) = — exp(—x) (x^ + 3x^ + 3x) 

48 

Inserting the results of Eq.(^) into Eq. (^), one can cast Tfi in the form: 


( 36 ) 


Tf^ = 


2G^ 


eff 


V2M, 


N 


47riin 


(37) 


where is : 


kl 


1 ( (^{k) • fki (J{i) ■ rkl _ (k ), (i) 


. 0-^0 


(38) 


kl / Z, Uij^ 

In Eq. (^) we have written M.fi in terms of the following known nuclear structure matrix 
elements W: 

Mgt,n = {{F\ ^T^'*T^^cr{k) ■ cr(/) (—) Fn{xa)\I)) 

kl 

Mf,n = ((Fl^rf rf (^) F^{xa)\I)) 

H ^^kiJ 

Mgt' = {{F\Y,Tf'^Tfaik) ■ (Til) (—) [^Fa{xa) - Fb{xa)] \I)) 

kl ^'^klJ 

Rn 


^F' = ((E| [3^4(xa) - Fb{xa)] |/)) 

\TklJ 

(fc) (0 ( ^jk) • ^kl Ojl) ■ Tkl _ (k) _ m 


kl 


vIt' = {{f|E4“'4‘'I4^II3 

kl ^ ^kl 


kl 


FBixA)\I)) (39) 


The numerical values for the case of ^®Ge quoted in ref. are given in Table I for easy 
reference. We can decompose Ali?/ as : 


Mfi = mPj + Mf} 


Note that in ref. a slightly different notation is used; E 4 = 3Ea — Fb and E 5 = Fb- 
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( 40 ) 


and finally obtain : 


<1 = 
Mf! = 


9a 


Ml 


l/\Mw 
2 ^ 


m. 


- Mgt' 
Mgt,n- 



mPj = +8.15 X 10"^ 
= -6.27 X 10“^ 
Mfi = -5.45 X 10-2. 


(41) 


Let us note that in Eq. the part of the nuclear operator which is independent of A 
coincides exactly 0 with the result of ref. (Neglecting the time ordering in the hadronic 
current is equivalent to the limit A ^ 0). Our result in Eq. (|37|) shows the expected 
scaling ~ 1/Mjy with the composite neutrino mass Mjy; exchanging a very heavy Majorana 
neutrino one expects a factor Mjy x from the neutrino propagator and the exchange 
of heavier particles reduces the probability of the decay. This behaviour does not appear in 
the formulae of ref. |^. 


V. DISCUSSION 


Calculation of the half-life of the decay from Eq. (|37|) is now straightforward. We could 
use the results of ref. 0 but prefer instead to use the already well known phase space 


factors given in which take into account the distortion of the electron wave functions 
going beyond the Rosen-Primakoff approximation. The inverse half-life (T 1/2 = log2T) is 
given by: 

1 f „ .. „ „ „ „ . d^Pi d?P2 


rp-l _ 

-‘- 1/2 — 


log 2 

From Eq. ( 0 ) one has: 


I \Tf,\^27iS{Ei-EF-Ei-E2) 


(27r)32Ei {2'Kf2E2 


(42) 


^apart from a normalisation factor. 
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IT’/iP = E \HP2){vlL + vlR)i’^{pi)[' 

^ 0 electron spins 

The electronic wave functions are given by: 

'lp{Pl,2) = Fo{Z + 2, Ei^ 2) u{pi,2) 


( 43 ) 


(44) 


where Fq{Z+2, -^ 1 , 2 ) is the well known Fermi function describing the distortion of the electron 
wave in the Coulomb field of the nucleus. We have then with straightforward algebra: 


E 

electron spins 


\^iP2){vlL + r|lR)^p^ = 


— Fq{Z + 2, El) Fq{Z + 2, E 2 ) 

electron spins 

= Fq{Z + 2 , El) Fq[Z + 2 , E2) {rj\^ + ? 7 ^) 2pi ■ p2 


E HP2){vIL + vIR)u^{pi)\‘‘ 


(45) 


As regards the phase-space integration, we adopt the standard notation and define: 
{Gf cos 6cgA)'^ml 


^Ou — 

Goi = 


647r^ 

•^Oi' 


2P1E1P2E2 


\og2 {meRof 


m'i 


S(^Ej — Ep — El — E 2 ) X 

Fq[Z 2, El) Fq[Z 2, E 2 ) dEi dE2 


(46) 


The phase-space integral in Eq. (^) is known in the literature and its value, which completely 
takes account of the Fermi functions, is Gqi = 6.4 x (see ref. [^). As for the 

half-life, we can write it as: 

(47) 

Experimentally, QulSf] decay has never been observed so far and therefore we have avail¬ 
able only lower bounds on the half-life bound. 

(48) 


Ti /o > Ti 


lower bound 


We can use the constraint in Eq. 
compositeness parameters: 


1/2 > -1-1/2 


into Eq. 


in order to obtain constraints on the 


/ 

< M'f 

( 

1/4 

lower bound i ^4 \ 

<^01 ^1/2 [Vl + Vr) 

Ac 

[m\) 

\Mfi 

1/2 


- 1/4 


( 49 ) 
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Let us now consider the neutrino-less double beta decay: 


76 


Ge 


,76 


Se -f- 2 e 


( 50 ) 


for which the HEIDELBERG-MOSCOW (3(5 experiment provides the current lower bound 
on the half-life P2i: 


> 7.4 X lO^V 

Substituting this value and the other numerical constants (t/l = rjR = 0) into Eq. 
obtain hnally: 

1,1 ,Qn„ V /jVv V'" 

- '®®lTeVUTeV/ ' 


( 61 ) 


we 


(52) 


Eqs. ^ ^ are the central results of this work. They describe the constraints on 

compositeness parameters imposed by the non observation of neutrino-less double beta decay. 
A few comments are in order. 

We have derived, according to the closure approximation, an expression for the nuclear 
matrix element of the form Aipj = + ■M.^pi where coincides with the matrix 

element given in ref. 00. and Mpj is the hrst order correction in powers of A. Thus in 
the limit A —0 we recover the previously published result []7|,p!2|, obtained neglecting the 
time ordering of the hadronic charged current c.f. Eq. p. 

Another important improvement of the present work with respect to is that the 

matrix element has now been expressed in terms of known matrix elements 0, so that 
while previously only a bound |A4^/| < 0.34 could be given, we have now a dehnite number 
for MFI c.f. Eq. Thus while in [^|T^ only a most stringent bound for the compositeness 

parameters could be quoted, here we have derived the bounds imposed by 0i'(3(3 non obser¬ 
vation. For dehniteness we have shown in Fig. ^the present bounds on |/| for (Ac = 1 TeV) 
as function of the heavy neutrino mass The bound (solid line) is compared with the 
most stringent bounds quoted in [^,|^]. We can see that the true bound is, as expected, 
somewhat weaker than the previously published one. The same conclusion can be drawn 
from Fig. where we give, as a function of Mpf, the lower bound on Ac for (|/| = 1). 
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The fact that within the closure approximation we obtain an analytic expression of the 
nuclear operator identical to the previous one (c.f. |]7|JT^) up to a term linear in A, which 
turns out to be numerically non negligible, gives us conhdence on the correctness of our 
result, namely Eq. (|i0|) and Eq. (^91) . 

We believe that the analysis given by Takasugi [Q, does not support his conclusion given 
in Eq. That result is in contradiction with the physical expectation that exchanging an 
heavier Majorana neutral reduces the probability of the neutrino-less double beta decay and 
thus must impose weaker bounds on the remaining parameters. We give, in the appendix, 
a detailed discussion comparing our calculation with that of Takasugi which show explicitly 
why he obtained the wrong scaling behaviour. 

A remark on the calculation of the matrix element is also due. While we think we 
have proved the correctness of the analytic expression for the nuclear operator involved 
in the Oi//3/3 with a composite Majorana neutrino, there still remain some uncertainty on 
the actual value of the matrix element Aipi- As reported in ref. |^, WIf' and M.gt' are 
quite sensitive to the value of m^, the cutoff parameter in the nucleon form factor, and Iq, 
the cutoff parameter of the short-range correlations between nucleons. This is due to the 
fact that the radial functions appearing in A4f' and M.gt' Eq.(|59|) are not positive dehnite 
resulting in a delicate cancellation within the radial integrals. This does not happen for 
M-t'-i ■M.gt,n and M.fn which are rather stable in the region 650 MeV < < 1.5 GeV. 

The numerical values used here to get Eq. (^) refer to the values = 0.85 GeV and 
Ig = 0.7 fm. 

The sensitivity of some of the nuclear matrix elements involved in this calculation to¬ 
gether with the fact that M.^fi is also quite important simply indicates that due to the 
coupling appearing in our effective Lagrangian, nuclear physics aspects of the neutrino-less 
double beta decay calculation become more important. Further investigation is necessary. 
In particular, we believe that the term I' in Eq. p8| ) containing in the numerator, gives a 
bigger weight to the high go region and in order to properly account for the nuclear physics 
it might be necessary to go beyond the closure approximation. 
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Let us now compare our result in Eq. (|5^) with those from high energy experiments. 
The ZEUS and HI collaborations (DESY) have recently published p3| , p^ results of a direct 
search of singly produced excited states in electron-proton collisions at HERA. They have 


studied the reaction ep —>• l*X with the subsequent decay I* —>• I'V where V = ■j, Z,W 
(see Fig. 6). Upper limits for the quantity \J\cvi*e\^ + \dvi*e \^/x Br^'^^(/* —> I'V) are 
derived as a function of the excited lepton mass and for the various decay channels. 
These experiments were sensitive up to 180 GeV for m^* and up to 250 GeV for rUe* {rriq*). 
For the purpose of comparing our analysis of double-beta decay bounds on compositeness 
with the high energy bounds, we quote here the limit on the v* coupling that the ZEUS 
collaboration has obtained at the highest accessible mass {my* = 180 GeV): 

\J\cWu*e\^ + \dwh 


^ u*e 


Ac 


X BrB2(iy* ^ tjW) < 5 X 10-2 GeV 


-1 


(53) 


Let US emphasise that these limits depend on the branching ratios of the decay channel 
chosen. For m^* = 180 GeV (the highest accessible mass at the HERA experiments p3|j2^ 
with Ac = 1 TeV, Br{y* —>• vW) = 0.61 and \cwu*e\ = \cwu*e\ one has: 


I/I <61. 


HERA 


(54) 


For the same values of m^* = Mjci and Ac one obtains from the Oi//?/? constraint i.e. Eq. 


< 3.84 




(55) 


We can thus conclude that the bounds that can be derived from the low-energy neutrino-less 
double beta decay are roughly of the same order of magnitude as those obtained from the 
direct search of excited states in high energy experiments. We also note that, in contrast 
to bounds from the direct search of excited particles, our QuPP constraints on Ac and |/| 
do not depend on any assumptions regarding the branching ratios of the decaying heavy 
particle. 

Finally, let us conclude by recalling that with respect to the work of ref. 00. we have 
improved the calculation of the phase-space using the exact values of the integrals given in 
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ref. 1^, thus removing the Rosen-Primakoff approximation for the Fermi functions. This 
however does not give appreciable changes in the numerical results as shown in Figures 2 
and 3. The results appear stable. 


APPENDIX: DETAILS OF CALCULATION 


Here we show in some detail where the calculation of ref. differs from ours and why 
the author of reference got the wrong scaling behavior for Eq. (2). The discrepancy is 
that the author there used throughout an effective four-fermion interaction, i.e. an effective 
Fermi theory. This is clear in view of the absence of W gauge boson propagators in his 
formulae. The use of an effective four-fermion interaction of course amounts to let Mw —> oo 
after having introduced the effective Fermi coupling constant Gp. Normally this is a good 
approximation if there is no other mass scale comparable with or larger than Mw, as is 
the case, for instance, in typical low energy processes. However, in our problem, another 
mass scale enters the game, namely the heavy neutrino mass which we assume to be 
much greater than Mw- It is then inconsistent to let Mw oo while keeping M^ hnite. 
Nevertheless, this was done in ref. p. Since we are to evaluate effects due to the heavy mass 
Mtv, we must include effects due to Mw- They can not be discarded a priori if we wish to 
discover the correct scaling behavior. 

Taksugi’s Eq. (7) in ref. [Q reads: 

exp [-iq ■ {x - y)] 


O ■ Glff 
% = ^—— 


(271) 


dxdy / dqM^ 


q^-M% 


{Nf I T(Jl(x)Jl(y) I N,) X 


^ (1 - P 12 ) P"(Pi, P 2 , q\ q) exp [i(E^x^ + E^y^) 


E^{E,, E2, q\ q) = ME2)a>^^aP^(ylL + (E,)(q^ - E,g^o)(q, + E^Qpo) 

f ' 


(Al) 


^eff - 


Af 


Gf cos 6c (homodoublet model) 


(A2) 


In Eqs. (A1-A2) Takasugi is also neglecting the momenta of the outgoing electrons ev¬ 
erywhere except in the electron wave function (which includes the relativistic Coulomb 
corrections of the nuclear held). Our Eq. (8) on the other hand may also be written as: 
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X 


Sfi — i {Gp cos 9c) 


2 l± 
Uc, 


(1 - -P 12 ) 


d'^x d'^y 


d^q exp [—iq ■ (x — y)] 


(27r)^ 


Ml 


K 


{Np I T{ji:{x)jay) I Ni) 


'4^{P2)(r^,x(Tup{vlL + X 


[(g - Pi)2 - M^] [(g + - M^] 

exp [i{pi ■x + p 2 -y)] (g - Pi)^(g + P 2 )^ 


(A3) 


In order to compare Eq. (A3) with Eq. (Al) we neglect, at this stage, in Eq. (A3), the 
electron momenta everywhere but in the wave function (as opposed to section III where 
we did so after having extracted the four-dimensional delta-function of energy-momentum 
conservation). Hence Eq. (A3) reads : 


5^. = i GIjjMn / d^x d^y 


exp [—iq ■ [X 


(27r) 


gZ 


Ml + ie 


y)] j^4 


{Np I T(j{:(x)jay) I Ni) 


(g 2 


Mlr + te)^ 


X 


(1 - P 12 ) V(Ei, E 2 , g°, q) exp i{Eix^ + Pag”) 


(A4) 


Our result coincides with Eq. (Al) in the limit Mw 00 . We have thus shown that our 
equations would coincide with those of ref. @] if we were to adopt the Fermi’s effective 
theory. 

Let us comment further on the calculation presented in ref. 0. From Eq. (Al) Takasugi 
hrst performs the g° integration picking up contributions only from the Majorana neutrino 
pole. However, for us, keeping the W propagator is essential. 

Let us now continue the calculation from Eq. (A4), but following the method employed 
in ref. [^, deviating somewhat from that of the present work, but with the exception of 
postponing the g° integration. Expanding the T-product, inserting a complete set of in¬ 
termediate states between the hadronic current operators, and integrating out xq, go, after 
some algebra (exchanging x <->■ g, g <->■ —g, /i i/ in the second term) Eq. (A4) yields : 


Sfi — 2 tt 6{Ej — Ep — El — E 2 ) Rfi 

Rf, = Ma, I d?x Xy 


d*q 


(1 - -P 12 ) 


(2nY 

tpuiq'i Eli E 2 


(qo + - ie) (g2 - Ml + ie){q^ - M^ + ie) 


exp [iq ■ {x - y)] {Np\J^{x)\X) {X\r{y)\Ni) x 

(AS) 
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where — Ex — Ep — Ei — Ex — ^j^iEi E Ep^ — 1/2(h/i — E 2 ') ~ Ex — \j2,(^Mj Mp'j. 
As in section III of this work, Ex is replaced by an average excitation energy {Ex) (closure 
approximation): 


^ A = {Ex) -^{Mi + Mp) ^ 10Mev. 


(A 6 ) 


Neglecting the electronic energies also in the numerator (in the tensor and using the 
spinor identity given in Eq. ( 0 ) we obtain : 


Rfi - 


72 


Mx J d?x(fy —^ exp [iq ■ {x - y)] {Np\J^^{x)r{y)\NI) x 


( 2 ' 


X H‘2)i.vlL + r/^i?)?/^(l) (A7) 


(go + A - ie)(g2 - + ie){q‘^ - + ief 

Now we make use of the non-relativistic impulse approximations for the hadronic charged 
current, just as in section III (inclusion of nuclear form factor etc.) and neglect terms in 
qoQi since we only want to consider O’*" —> O’*" transitions. We End: 

Rf,{0+^0+) = ^^^ij{2){r]lL + r]lR)^^{l)MNM^ [ 7 ^^ /l(g^) x 


72 

^ {Np\ exp {iq ■ rfc;)r|^^ q^I{q^ 

k,l ^ 

2\ Rk) Jl) 


(27r) 


E9a 


-V^^q^ - qW X{q^) a)' 


, gla^^^-a^^H'{q^)\\Ni) (A 8 ) 


where we have dehned: 


= 


= 

J'{q^) = 


d 




J{e 


f dqo 

\ 1 1 

1 

J 271 

.go - A + ie. 

{ql -uj% + ie) {ql - + ie) 


d 


J'{q^^ 




dqo 

271 


ql 


{ql -ujI + ie) {ql - + ie) 


.go - 2 ^ + R- 

with uix and ujw are the same as in section III. 

Eq. (A 8 ) is found to be equivalent to Eq. (^9]) because of the identity: 


(A9) 


dqo 

271 


go + A - ie 


2 ^ 1 /■ dqo 


2iA 


ql - IS? + ie 


f{ql) 


(AlO) 
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which implies that the integrals X, X' are proportional to I, I' dehned in Eq. section IV, 
namely : X = {i/2) I; X' = {i/2) I'. The fact Eq. (A8) coincides with Eq. (p9|) concludes our 
proof that including the W boson propagators in the calculation of ref. gives the correct 
scaling behaviour as found in the present work (and also found in ref. |]^). 
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TABLES 


TABLE I. Nuclear structure matrix elements for ^®Ge as reported by the authors of ref. [20|. 


Mgt,n 

Mfn 

A4gt' 

A4 F' 

A4f' 

1.13x10"^ 

-4.07x10-2 

-7.70x10-3 

3.06x10-3 

-3.09x10-3 
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FIGURES 


FIG. 1. Neutrinoless double beta decay (AL = +2 process) mediated by a heavy composite 
Majorana neutrino. 


FIG. 2. New bounds on the parameter |/| (for Ac = ITeV) from Qh'f3f3 (solid line) compared 
with the estimate, based on a upper bound of the nuclear matrix element {\Mfi\ < 0.34), given 
in ref. [p|,p^ (dashed line). The dotted line is the calculation of ref. augmented by the exact 

phase space calculation (no Rosen-Primakoff approximation for the Fermi functions). 


FIG. 3. New bounds on the parameter Ac (for |/|=1) from OuPP (solid line) compared with the 
estimate given in ref. |@JT^ (dashed line). The dotted line is the calculation of ref. |@jl^ augmented 
by the exact phase space calculation (no Rosen-Primakoff approximation for the Fermi functions). 
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Upper Bound on |/| 


Figure 2 
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